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In this work, we calculate the exact asymptotic quantum correlations between two interacting 
non-resonant harmonic oscillators in a common Ohmic bath. We derive analytical formulas for the 
covariances, fully describing any Gaussian stationary state of the system, and use them to study 
discord and entanglement in the strong and weak dissipation regimes. We discuss the rich structure 
of the discord of the stationary separable states arising in the strong dissipation regime. Also under 
strong dissipation, when the modes are not mechanically coupled, these may entangle only through 
their interaction with the common environment. Interestingly enough, this stationary entanglement 
is only present within & finite band of frequencies and increases with the dissipation rate. In addition, 
robust entanglement between detuned oscillators is observed at low temperature. 
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Over the last decades, non-classical correlations such 
as discord and entanglement have been widely acknowl- 
edged to play a central role in quantum mechanics. In 
particular, the robustness of entanglement against de- 
coherence in continuous- variable (CV) systems has been 
object of study in a number of recent publications [TUTU]. 
Due to the fact that entanglement is a valuable resource 
in quantum communication and information processing 
with CV, the investigation of the conditions under which 
bipartite states indefinitely preserve their non-separable 
nature in presence of noise and dissipation is not only of 
fundamental importance, but also of practical interest. 

Gaussian two-mode states occupy a priviledged posi- 
tion among all entangled CV preparations, since they 
combine easy experimental realization (e.g. the output 
beams in an optical parametric oscillator), with simple 
mathematical description in terms of their second order 
moments, which are experimentally determined via ho- 
modyne detection [11 j. These CV systems also find phys- 
ical realization in mechanical oscillators, such as trapped 
ions coupled by their electrostatic interaction [12]. A 
simple but rather general model describing those systems 
consists of two coupled detuned harmonic oscillators, lin- 
early interacting, and in contact with a bath. 

The stationary entanglement of two resonant modes 
in a squeezed state under decoherence was studied in 
Refs. [2 , whereas dynamical features of the non-resonant 
case were addressed in [3 . Additionally, different long- 
time behaviours were experimentally observed in [13] af- 
ter simulating decoherence on one of the modes. Other 
recent publications such as [3 [8] treat the dynamics of 
quantum correlations in similar systems for a variety of 
structured environments. Nonetheless, it must be noted 
that the majority of these works rely on weak system- 
environment coupling assumption. 

On reference to quantum discord [14 , it was intro- 
duced to capture the quantumness of correlations. Since 



the majority of quantum states, including most of the 
separable ones, have non-zero discord [15], it is said to 
be more general than entanglement as a measure of non- 
classicality. Discord has recently attracted attention in 
the field of quantum information, now that for some com- 
putational models wich make no use of entanglement, it 
is believed to provide exponential speedup over their best 
known classical counterparts [T6] . 

In this Letter, we provide exact analytical formu- 
las fully characterizing the Gaussian stationary state of 
two interacting non-resonant modes in contact with a 
common thermal environment by solving the quantum 
Langevin equations in the asymptotic limit [4] [5]. We 
then compute quantum discord and entanglement ex- 
actly, being our main interest the robustness of stationary 
quantum correlations against temperature and dissipa- 
tion, in the less studied case of finitely detuned modes. 

Our system consists of two harmonic oscillators 
coupled through a term bilinear in the coordinates 

H s = Ei=i f + h"i x i + I 0*1 " x *) 2 ' The in- 
action with a thermal environment can be modeled 

by coupling them to a common bosonic heat bath, 
which yields: H = H s + £ M + + 

T Oi, x 2 ) E/i^/i + J r (x 1 ,x 2 ) 2 ^ 2 /2 [17, where Q 2 is 
given by 2 duo J (uj) /oj : and J (uo) is the bath spectral 
density [18] . We use units in which h = ks = 1 • 

As noted in Refs. [2j [7] among others, under suitable 
conditions, an initially separable state may become en- 
tangled in this model, even if k = 0. However, this ability 
of the common environment to create non-classical cor- 
relations between the two parties highly depends on the 
physical distance r that separates them [4[ [10] , so that 
distant dissipative oscillators would be better modeled 
by considering independent identical environments. 

In the following, all averages are taken with respect to a 
separable initial state in the system's and bath's degrees 
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of freedom p (0) (g) pr where p (0) is a Gaussian state of 
the two oscillators, and pT is a thermal equilibrium state 
of the environment at temperature T. 

By choosing T x 2 ) = #1 + ^2 in # one may derive 
the following quantum Langevin equation [T8] : 

2 

+ (cj? + ^ 2 ) Xi + ^ [fc (xi - Xj) + ^ 2 x J -] 5^- (1) 
i=i 

2 ft 

i=l J -°° 

where Sij = 1 — ^ . The quantum force F (t) acting on 
both oscillators and the dissipative kernel x (t) are con- 
nected by the Kubo relation x (t) = i0 (t) ([F (t) , F (0)]), 
where 6 (t) stands for the Heaviside step function. For 
this choice of F(xi^x 2 ), the overall system Hamiltonian 
H is quadratic in positions and momenta, so that for 
uji 7^ co>2 the asymptotic state (which is independent of 
p(0)) is guaranteed to be Gaussian. In the case of res- 
onant oscillators, the normal mode decomposition leads 
to an effective decoupling of the relative degree of free- 
dom from the bath, and thus the stationary state would 
depend on the initial condition [2 . In the following we 
shall restrict ourselves to non-resonant oscillators. 

We can solve Eq. by Fourier transforming and 
denoting f (uj) = J^dt e iu)t f(t). Due to the lin- 
earity of the problem, this formally yields: X\ (uj) = 
C6ij (uj) F (uj) , where the generalized susceptibilities 
ctij (uj) are elements of the inverse of the matrix: 
(a -1 )^ = uj 2 + Q 2 — uj 2 + k — x(uj) and (a~ x )ij = 
— k + ^ 2 — x(uj) for i ^ j. 

The central object of our study is the asymptotic co- 
variance matrix T, defined as = (R%Rj) — [Ri, Rj] /2, 
with R = {x±,pi, ^2,^2} which, up to local displace- 
ments, fully characterizes any Gaussian state of our bi- 
partite system. Provided with the susceptibility matrix 
we may compute the power spectrum (xiXj) = 

dt e iu; ( t_£/ ) (xi (t - t') Xj (0)) of the stationary two- 
time correlation (xi (t — t') Xj (0)) = (xi (t) Xj (t f )) in 
terms of the power spectrum of the environmental force: 



( X i X j)«, = J2 k 1 aik ( W ) ( FF ^ 



(2) 



The rest of power spectra (RiRj)^, can be readily ob- 
tained from Eq. pi) since (piXj) 
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-iuj (xiXj) and 

(PiPj) = ^ 2 ( x i x jT • The full asymptotic covariance 
matrix then reads Tij(t,t) = f^° oo (duj/27r) (RiRj) u — 
[Ri, Rj] /2. 

To proceed further in deriving analytical expressions, 
we shall choose an Ohmic spectral density with Lorentz- 
Drude cutoff J (uj) = 2juj/tt(1 + uu 2 /uj 2 ), that leads to 
a dissipative kernel x ( u ) = 2juo 2 /(uj c — ico), and a bath 
correlation function (FF) u = ttJ (uj) [1 + coth (u;/2T)], 
where uj c is the cutoff frequency and 7 the dissipation 
rate. The renormalization constant appearing in Eq. ([I]) 



is Vt 2 = 2juj c for this choice of J (uj). We start from Eq. 
([2} to get 

where gij (uj) and h (uo) are 5th order polynomials de- 
pending on cji,cj2,^ c ?7 and k. After a lengthly calcula- 
tion, these integrals can be solved to yield [19]: 

( Xi (t) xj (t)) = (xiXj) T ^ OQ - a k ,j) ^ ( 1 - ^ ) 

k=l ^ ' 

^ (t) Pj (t)> = (p^)^ - J2 (1 - i^f) > 



k=l 



while the remaining convariances are all zero. cl£^ and 
bfr'^ are temperature-independent coefficients for which 
we have explicit formulas, z k are the five complex roots 
of h(uj), and ijj stands for the logarithmic derivative of 
Euler's gamma function. (RiRj) T ^ oo denotes classical 
mean values, to which the summations are quantum cor- 
rections. 

To gain some insight into the physics of the prob- 
lem, we will set k = and place ourselves in the 
limit of 7/w c <C 1. To first order in this small pa- 



rameter, the roots z k become zbcoj ± 
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2uj 2 ) + Additionally, in 



the very low temperature regime (i.e. T <C uJi) the 
asymptotic expansion ip (1 + z) ~ \ogz + (2z)~ may be 
justified. Inserting all this into the formulas of the T^-, 
and retaining terms only to first order in T /uj Ci uJi/uo c and 
j/uo c one gets the temperature-independent expressions: 
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By setting 7 = 0, the variances (x 2 ) and (p 2 ) become 
those of the ground state of the corresponding oscillators 
without dissipation, and the asymptotic (separable) state 
has a diagonal covariance matrix. Taking into account 
that uJc/uJi ^> 1, one sees that increasing 7 leads to a 
reduction of (x 2 ) while (p 2 ) gets larger, as expected. 
The terms proportional to log uo c /uJi in (PiPj), also occur 
in the same region of parameters for a single damped 
oscillator [T8] . 

These analytical and exact expressions for the asymp- 
totic covariances T^- are one of the main results of this 
Letter. From them, we shall compute all quantum cor- 
relations without any assumptions on the different time 
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FIG. 1. (a) Asymptotic entanglement and (b) discord in the 
weak coupling regime as functions of k and UJ2- We have 
chosen cji = 10, 7 = 0.01, uj c 500 and T 0.5. Both 
measures of quantum correlations exhibit a similar behavior: 
They increase with 002 /k — >• 0. Dashed lines are contours and 
the color legends are plotted in the upper left corners. 

scales involved in the problem. The logarithmic negativ- 
ity [20] is a suitable measure of entanglement is this case. 
It is based on the positivity of the partial transpose sep- 
arability criterion [21 , which is necessary and sufficient 
for bipartite Gaussian states [22] . When it comes to dis- 
cord, it is esentially the difference between the mutual 
information of the bipartite state, before and after per- 
forming a complete measurement on one of the parts (in 
what follows, mode two). Even though its calculation is 
usually very involved, if one restricts to Gaussian gener- 
alized measurements it is straightforward to compute for 
bipartite Gaussian states, as recently shown in [8]l23ll24]. 
We now turn to present the results thus obtained: 

Weak dissipation. — We will look first into the weak 
dissipation regime, where the relaxation time scale tr = 
7 _1 is much larger than the bath correlation time, given 
by tb = V\ X whenever uj c > v\ [25]. We shall also pick 
a low temperature T for which non-zero asymptotic en- 
tanglement might exist. Results are presented in Fig. [I] 
In accordance with [2] , we observe Ej^ ^ out of reso- 
nance at low temperatures, and overall similar behaviour 
of entanglement and discord. 

In the limit k ^> u^, the oscillators, whose effective 
frequency Co 2 is given by uo 2 + k + 2~fuo c1 become nearly 
resonant and thus, approximately decouple in the vari- 
ables r]± = (xi d= X2) /y/2 and tt± = (pi ±^2) /V%- As 
ujk -+ 0, {( V 2 + ) , (tt 2 _ >} ->■ 00 while {(n 2 + > , {if_ )} -> 
[5]. This corresponds to a non-symmetric two- mode 
squeezed thermal state with infinite squeezing r, i.e. the 
ideal EPR (maximally entangled) state [26] . 

Strong dissipation. — We achieve this regime by in- 
creasing 7, so that the relaxation time tr becomes com- 
parable to the bath correlation time tb- Results are pre- 
sented and summarized in Fig. [2] 

Even at k = 0, there may exist non- vanishing station- 
ary entanglement. The corresponding asymptotic states 
are thus non-Gibbsian, as could be expected at suffi- 



FIG. 2. (a) Asymptotic entanglement in the strong dissipa- 
tion regime as a function of k and CJ2- Due to the bigger 
decoherence entanglement is mostly lost, though its maxi- 
mum is still attained as U02 /k — »• 0. Interestingly, in the 
inset of (a), we can observe a small region around k — 
where non-classical correlations are enhanced by the dissipa- 
tion. This purely non-Markovian effect, enters into our equa- 
tions through the renormalization term Q 2 xiX2 in the Hamil- 
tonian H, which for sufficiently high dissipation rate, results 
in an effective environment-assisted coupling between the os- 
cillators, (b) Asymptotic entanglement (continuous line) and 
discord (dashed line) as functions of UJ2 alone for k = 0. For 
both plots, parameters are the same as in Fig. [I] except for 
the dissipation 7 = 0.5. 



ciently low temperatures [27 . Actually, one sees in Fig. 
[2] (b) that for fixed cji, 7, T and co> c , there exists a band 
of frequencies 002 £ [uj',uj"] for which the asymptotic 
state is non- separable. As the temperature increases, the 
amount of entanglement gets reduced, while uo' and uj" 
come closer to each other. On the contrary, by increas- 
ing the dissipation rate 7, the maximum attainable Ej^ 
and the bandwidth uj" — uo' get larger. Reducing the fre- 
quency uji also leads to narrower bandwidth. One may 
also see that uj' is almost insensitive to changes in o; c , 
while uo" slightly decreases as uu c grows. 

The simple model H s = EM/ 2 + K 2 + ^ 2 K7 2 ] + 
Q 2 xiX2 proves useful to understand these features. Note 
that the renormalization Q 2 = 27Co> c has been intro- 
duced to account for the bath-mediated coupling. If 
Wi/Q 2 <C 1, the oscillators are again near resonance 
and approximately decouple in the coordinates {rj±, tt±}. 
Then, E M ^ whenever v 2 _ = <77?_> <7t^_) < 1/4 

It is easy to see that the lowest symplectic eigenvalue of 
the partially transposed covariance matrix z>_ grows with 
T and decreases with Q 2 and coj. Thus, an increase in the 
temperature must be compensated by an increase in uj' 
for the inequality to be satisfied. Conversely, an increase 
in either Q 2 or uj\ allows for a reduction of the critical 
frequency uj' . The existence of an upper bound uo" for 
this non-separable region is related to the structure of 
the bath, so that its behavior can not be inferred from 
this simple non-dissipative model. 

Finally, we also remark that in our symplified model, 
V- decreases as \Q 2 — k\ grows, since the effective coupling 
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discuss the non trivial structure that quantum discord 
shows when the stationary states are separable. 
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FIG. 3. The asymptotic discord displays a very different be- 
haviour in the strong dissipation regime. Since the entan- 
glement has considerably decreased, it no longer dominates 
over other quantum correlations and, as could be expected, 
remains non-zero even when the asymptotic state is separa- 
ble. It displays a maximum around k c± 20, marked with the 
dotted line, and also vanishes in the limit {/c,^} — > (see 
text for discussion). All parameters are the same as in Fig. 

ca 



between the oscillators is the net result of two compet- 
ing effects, that may eventually cancel each other: The 
environment- assisted coupling and the direct interaction 
between the oscillators. This explains qualitatively why 
E^f = for intermediate couplings in Fig. |2j 

Finally, Fig. [3] contains results on discord in this 
regime. The global purity \i may be found to be roughly 
constant for all stationary states, except those with 
{/c,^} — » 0, while z>_, is maximum on the dotted line 
of figure [3| This suggest that at fixed /i, Gaussian quan- 
tum discord increases with z>_, which is consistent with 
[23] , where this was seen to hold for squeezed thermal 
states. Additionally, since the partial purity ji2 — >• as 
{/c,^} —> the marginal stationary states p2 = Tri(pi 2 ) 
become maximally mixed in this limit, therefore yield- 
ing a quantum-classical state [28] with zero discord as 
revealed by measurements on mode two. 

At higher temperatures and stronger dissipation rates, 
the fate of the asymptotic entanglement is to disappear 
completely, in contrast with the asymptotic quantum dis- 
cord, which is always non-zero as could be expected [T5] . 

To summarize and conclude, we have exactly solved 
the problem of the asymptotic quantum correlations in a 
system of two coupled non-resonant harmonic oscillators 
in a common bath. We provided exact analytical formu- 
las for the covariances, fully characterizing any Gaussian 
stationary state of the system, and used them to compute 
the asymptotic entanglement and discord, in the weak 
and strong dissipation regimes. We found that stationary 
entanglement between non-resonant oscillators may exist 
at sufficiently low temperatures. Furthermore, for a finite 
band of frequencies, it may be created only due to the 
environment- assited interaction between the two modes 
in the strong dissipation regime. Also in this regime we 
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